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ABSTRACT 



Several synoptic problems, previously solved by graphical overlays, can 
be solved with a programmable calculator. The HP-41CV hand-held calculator is 



used to determine surface winds, the 
regional area (Monterey, California), 
condensation level from surface data. 



probability of average rainfall for a 
contrail formation, and the lifting 
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I. Introduction 

This report describes four programs that use the HP-41CV calculator’s 
programming capability in dealing with interesting meteorological problems, 
including determining surface winds, the probability of average rainfall for a 
regional area (Monterey, California), contrail formation and the lifting 
condensation level from surface data. 

The programs were prepared to use the hand-held calculator in solving 
synoptic problems. The first problem is the estimation of surface wind at a 
point from a regional sea-level pressure analysis or forecast. This problem 
arose in the synoptic laboratory where wind forecasts were prepared with 
numerical weather prediction guidance (facsimile sea-level pressure maps). 
The immediate solution is a geostrophic wind scale applied to the facsimile 
chart. The HP program provides an alternative solution. The program 
calculates the geostrophic and gradient wind, as well as the estimated surface 
wind, using Ekman boundary layer assumptions, by simply entering latitudes, 
longitudes and pressure readings from a sea-level pressure chart. 

The second problem involves preparing a statistical precipitation 
forecast for a set of graphical tables. Using a tutorial program in a hand- 
held calculator, the need for the graphical tables is eliminated. The program 
determines the probability of average rainfall over the Monterey, California 
area for the period 0800 PST forecast day to 0800 PST of the following day, 
for the wet-season months, November through April. This procedure is based on 
the Renard (1972) graphical objective technique for forecasting 24-h rainfall 
at Monterey, California. From the data entries of 500 mb heights used to 
calculate geostrophic relative vorticity at Monterey and a point 8° latitude 
upstream from Monterey, as well as sea-level pressure data, the program uses 
a series of pre-established data tables to arrive at a probability of 24-h 
rainfall for Monterey, California. 
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The third problem involves forecasting contrails from upper-level 
sounding data. It was desired to have the HP make the decision from input 
sounding data as an alternative to plotting the entire sounding and using a 
graphical overlay. This program determines contrail formation at a given 
pressure level based on the curves found in the Contrails Forecasting Manual 
of Chief of Naval Operations (1964). A least squares approximation was made to 
the curves involving an expression for critical temperature which is a 
function of temperature, pressure and relative humidity. If the relative 
humidity is estimated (dew-point temperature is unknown), there is a series of 
criteria to determine whether or not contrails will occur. 

The final problem is to determine the lifting condensation level (LCL) 
without using a thermodynamic diagram. Inputing values of surface 
temperature, surface pressure and specific humidity to the HP results in a 
lifting condensation level. The LCL is the height at which a parcel, when 
lifted dry-adiabatically, becomes saturated. The program uses the Clausius- 
Clapeyron equation to arrive at a final expression for the LCL that is 
determined by values of surface temperature, surface pressure and specific 
humidity. 
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II. Computation of Surface Winds 

A. Objectives 

This program will give the forecaster, either ship-board or at a regional 
center, an easy and accurate method of computing the geostrophic, gradient 
and estimated surface wind from a MSL pressure analysis or model .forecast. 
This method is more complete than a graphical geostrophic wind scale overlay. 

B. Principles 

The geostrophic wind speed and direction are obtained from the east-west 
(ug) and north-south (vg) components. 



_1 3P 
fp 3y 





tan 9 e = — 

Ug 



where ug and vg are geostrophic wind components, 

f : 2 H sin <t> ; is latitude of forecast point, 



p = 1.3 kg/m 3 ; typical sea level density value 



3x, 3y = 1° lat = 111 km 



Geostrophic wind fVg ) computed by 
the finite differencing scheme 
with 1° latitude grid spacing. 



Figure 1 
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The user enters the latitude and longitude of the point (Northern 



Hemisphere) where a wind is desired. Next, the user enters four pressures 
that are 1° latitude north, south, east and west of the center point used in 
the finite differencing technique: 



The gradient wind speed is obtained from the geostrophic wind speed and 
the radius of curvature (R) . 



The user enters a radius of curvature (in kilometers) which is representative 
of the isobaric curvature in the forecast area. (See Figure 1.) 

In a cyclonic situation (R > 0), the magnitude of the gradient wind is 
less than the magnitude of the geostrophic wind (subgeostrophic) ; if R < 0 
(anticyclonic) , the gradient wind speed is greater than the geostrophic 
(supergeostrophic) . If the curvature is zero (R = ® ) the magnitude of the 
gradient and geostrophic winds is the same. Regardless of the curvature, the 
gradient wind is from the same direction as the geostrophic wind. 

The estimated surface wind speed and direction are obtained from the 
gradient speed and direction using Ekman layer assumptions. Over the land, 
the gradient wind speed is multiplied by a factor of 0.6 which is an 
approximate amount that the speed is reduced by friction. The resulting speed 
is an estimated surface wind speed. The actual direction is rotated counter- 
clockwise 20° from the original direction. Over the sea, the gradient wind 
speed is multiplied by .81 and direction is rotated by 10°. The oceanic and 
land frictional effects follow the observational summary of Sheppard (1969) 
and the discussion of Petterssen (1956) and Haltiner wand Martin (1957). 
Sheppard (1969) and Mendenhall (1967) indicate 10° is a reasonable estimate of 




(See Figure 1 . ) 



3y 222 km 1 3x 222 km 




where R > 0 cyclonic 

R < 0 anticyclonic 
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frictional veering over the ocean from weather ship data. Over land, the 
frictional effects are naturally heavily dependent upon surface roughness so 
the 20° veering is only a general estimate. 



C. User’s Guide 

1) Make sure the calculator and printer (NOR mode) are on. 



2 ) 



XEQ 



alpha I SIZE alpha 



will prompt SIZE 

3) Type 100. 

4) Enter magnetic cards as labelled 1-9 - - 
After the program (WIND) is "WORKING" and entered, it is ready to be 
executed. 

WIND 



5) 



6) 

7) 

8 ) 



XEQ 



alpha 



alpha 



The prompt CTR LAT? will appear. 



Press the R/S button to continue. 



(Do this each time you answer a prompt.) 

9) The prompt CTR LON? will appear. 

Enter the longitude of the center forecast point. 



10) Press R/S 



11) Answer the prompts for the four pressures (mb). 

(Remember to do this each time.) 



12) Press R/S 



13) The geostrophic wind speed (m/s) and direction will be displayed. 

14) Enter a radius of curvature (m) after prompt RADIUS? 

15) The gradient wind speed will be displayed. 

16) The prompt LAND or SEA? will appear. 

If the forecast region is over land, enter "LAND"; otherwise, enter 
"SEA". 
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17) The estimated surface wind speed and 



direction will be displayed and 



the program is done. 

18) To run the program again: 



XEQ 



alpha 1 WIND 



alpha 
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D. Examples 

To illustrate the use of the program the surface wind will be estimated 
for points A and B in Figure 2. The observed surface wind speed and direction 
for forecast point A at 0000Z on 26 November 1981 is 10 m/s from 315° . For 
point B, the surface wind speed is 20 m/s from 295°. 

Execution of the "WIND" program on the HP-41CV calculator results in a 
surface wind speed of 14 m/s from a northwest direction of 319° for forecast 
point A. Entering the latitude/longitude position (43N 130 W) along with four 
pressures 1° latitude north, south, east and west (1007 mb, 1012 mb, 1009 mb, 
1012mb respectively) results in a geostrophic wind of 20 m/s from 329° . A 
radius of curvature of 1,100 km allows a gradient wind to be computed (18 
m/s). Finally, since forecast point A is an ocean location, the oceanic 
Ekman layer boundary assumptions are applied. The resulting estimated surface 
wind is 14 m/s from 319° . This estimated wind is in agreement with the 
observed real wind speed (10 m/s) and direction (315°). 

For forecast point B (46N 135W) the four pressures (1008 mb, 1011 mb, 
1006.5 mb and 1013 mb) result' in a geostrophic wind of 24 m/s from 295°. 
From the radius of curvature (1,300 km), a gradient wind of 21 m./s is 
calculated. The final estimated wind speed for this ocean location is 17 m/s 
from 285° which compares with the observed real wind speed (20 m/s) and 
direction (295°). 
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Figure 2 

Surface conditions observed at 0000Z 26 November, 1981 for 

illustration of surface wind program. 
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Program 
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Surface Winds Program. 
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E. Program 







Table I 






Surface 


Winds Program. 




■ ♦LBL “HIND" 


/ 


101 PT&n 


151 8 TO 18 


FIX 0 


1.3 


102 188 


152 GTO -GRPIi- 


"CTp! LAT? 1 ' 


. / 


103 + 


153*LBL -EOST- 


Ha FKl.lHPT 


54 RCL 08 


104 STO 18 


154 90.O 


95 $70 8! 


55 / 


195 GTO -GPW 


155 STO 18 


06 -CTP LON’" 


56 STO 10 


106*L8L -WSH- 


156 GTO “GRAD - 


87 FPOrtPT 


57 RCL 09 


107 8 


157*l8L "GRAD” 


88 STO 82 


58 EHTERf 


188 enter: 


158 -GEO. BIR FROM 


09 -PPS HO. LOT'" 


59 xt2 


189 RCL 18 


159 PVIEM 


18 PROMPT 


60 RCL 10 


118 X>Y’ 


168 PSE 


11 STO 03 


6i m 


111 GTO “TRE“ 


161 RCL 18 


12 -PRS SO. LA T? ’ 


62 * 


112 1/X 


162 VIEW 18 


13 PROMPT 


63 SORT 


113 RCL 09 


163 PSE 


1“ STO 84 


64 STO 11 


114 * 


164 -RADIUS’" 


15 “PRS EAST LQT>- 


■ 65 “GEO. SPEEIi: ‘ jf5 ATAN 


165 PROMPT 


16 PROMPT 


66 OVIEW 


116 CHS 


166 STO 25 


17 STO 05 


6,7 PSE 


117 278 


167 8 


18 -PRS WEST LOT? 1 


■ 68 RCL ii 


113 * 


168 ENTER: 


19 PROMPT 


69 V I EM 11 


119 STO 18 


169 RCL 25 


28 STO 06 


76 PSE 


128 GTO -GRfiD- 


170 X<=Y’ 


21 RCl 03 


71 8 


121 ♦LBL -IRE- 


171 GTO -ANTI- 


22 ENTERt 


72 ENTER: 


122 1/X 


172 RCL 11 


23 RCL 04 


73 RCL 09 


123 RCL 09 


173 RCL 25 


24 - 


74 X=Y? 


124 * 


174 / 


25 108 


75 GTO -CHK- 


125 ATAN 


175 RCL 03 


26 * 


76 0 


126 CHS 


176 / 


27 2.22 E05 


77 ENTEPt 


127 270 


177 .25 


28 / . 


78 RCL 10 


128 * 


178 ♦ 


29 1.3 


79 X=Y? 


129 STO 18 


179 SORT 


38 / 


88 GTO -CHKS- 


130 GTO "GRAD" 


180 .5 


31 CHS 


81 0 


131 ♦LBL -CHK- 


181 + 


32 STO 07 


32 enter: 


132 8 


182 STO 26 


33 RCL 0! 


83 RCL 89 


133 ENTER: 


183 RCL 11 


34 ENTER? 


84 X>Y’ 


134 RCl 18 


134 RCL 26 


35 -SIN 


85 GTO -WEN- 


135 X<-Y7 


185 / 


36 1 . 4584 E-04 


36 0 


136 G?0 -NOR- 


186 STO 27 


37 * 


37 ENTER: 


137 180.0 


137 GTO "ACT- 


38 STO 88 


88 RCL 18 


133 STO 18 


188*LSL -ANTI- 


39 RCL 07 


39 I i)V 


139 GTO -GR0D- 


189 RCL 11 


40 ENTER? 


98 GTO -TMO- 


148+LBL -NOR- 


198 RCL 25 


41 RCL 08 


91 1/X 


141 368.0 


191 / 


42 / 


92 PCI 89 


142 STO 18 


192 RCL 08 


42 STO 09 


93 * 


143 GTO “GRQTi- 


193 / 


44 RCL 05 


94 ATAN 


144»LBL -CHKS- 


194 STO 29 


45 ENTEPT 


93 STO 18 


145 0 


195 RCL 28 


46 RCL 06 


96 GTO -GRfilr 


146 ENTER: 


196 ENTER? 


47 - 


97?L3L -TUO- 


147 RCL 09 


197 .25 


48 100 


98 1/X 


148 X<=Y> 


198 + 


49 * 


99 RCL 09 


149 GTO “EAST- 


199 STO 28 


50 2.22 E05 


100 * 


150 278.0 


200 0 



9 



III. Determining the Probability of Precipitation for Monterey, California 

A. Objective 

To determine the probability of 24-h rainfall for a regional area 
(Monterey, California) for the period 0800 PST forecast day to 0800 PST of 
the following day, for the wet-season months, November through April. This 
procedure is based on Renard’s (1972) honogram method and statistics. 

B. Principles 

General Forecast Procedure . The probability of 24-h rainfall is 
determined by values of geostrophic relative vorticity and sea-level pressures 
at different locations. 

A first forecast parameter is obtained by taking the sea-level pressure 
at Eureka, CA (PEureka) and the difference between the 500 mb 
geostrophic relative vorticity at Monterey ( Cm) and a point 8° latitude 
upstream from Monterey ( Cs). Enter these values into a pre-established 
statistical table, which is stored on the calculator, to arrive at a forecast 
variable, Yi . (See Figure 4.) 

A second forecast parameter, Y 2 , is the value that results from taking 
the sea-level pressure difference between Monterey and Eureka (APw-Eureka) 
and Monterey and Las Vegas (APm-lv). (See Figure 5). These two parameters, 
Yi and Y 2 , are entered in a third table (Figure 6) which results in a final 
value that corresponds to a probability of average rainfall for the Monterey 
Peninsula. 

500 mb Geostrophic Relative Vorticity . To calculate the 500 mb 
geostrophic relative vorticity, five 500 mb heights are needed to use finite 
differencing technique to solve (£m = g/f V 2 Z) where m is the geostrophic 
relative vorticity at Monterey, V 2 Z is the Laplacian of *the 500 mb height. 
For example, in Figure 3 at Monterey (Zm = 558 dm) the four 500 mb heights 
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are 552 dm (Zn : 500 mb height to the north), 566 dm (Zb), 560 dm ( Z<? ) , and 
563 dm (Zw) with a grid spacing of 6° latitude. This results in a geostrophic 
relative vorticity of 1.8 x 10 -5 sec' 1 . 



Figure 3 

Computation of geostrophic relative 
vorticity for Monterey (£ m) and a 
point 8° latitude upstream (£s). 



The 500 mb geostrophic relative vorticity is computed for a point 8° 
latitude upstream (8° corresponds to a typical daily progression of short- 
waves in the westerlies of 22 kt during the winter season) from Monterey in 
the same way. The difference between the two vorticities (£s° -£m ) is used 
as a forecast factor. 

Sea- Level Pressure . Yi , a first forecast parameter, is obtained from 
the computed vorticity difference and a sea-level pressure at Eureka, 
California (Figure 4). Y 2 , a second forecast parameter, is obtained from 
the difference between the sea-level pressure at Monterey minus the sea— level 
pressure at Eureka (APm-eui- e ka ) and the sea-level pressure at Monterey minus 
the sea-level pressure at Las Vegas (APm-lv) (Figure 5). 

Probability of Rainfall . From the two forecast parameters, Yi and Y2 , a 
value Y3 is obtained which corresponds to a probability of rainfall for the 
Monterey, CA area for the period 0800 PST forecast day to 0800 PST of the 
following day, for the wet-season months, November through April. (Figure 6.) 
The corresponding values between Y3 and a probability of average rainfall are 
on data cards that are entered during the program execution. 
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Figure 4 

Yi = fCPEureka, ?8 -Cm): first forecast 
factor of the raincaster. 
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Figure 5 

Ys = f ( APm- £u r i» k a , APm-lvj: 
second forecast factor 1 
the raincaster. 



Figure 6 

Y^ = f(Yi, Y?. > : third 

forecast factor of the 
raincaster. 
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c. 

1 ) 

2 ) 

3) 

4) 

5) 

6 ) 

7) 

8 ) 
9) 



User’s Guide 

Make sure the printer is attached and in the NOR mode. Both printer and 
calculator should be on and the calculator’s memory cleared. 



XEQ 



alpha SIZE alpha 



will prompt SIZE 

Type 035. 

Load the program cards (1-18). Continue after the prompt "WORKING." 



XEQ | alpha | RAIN alpha 



Enter 500 mb heights (dm) for each prompt. (For example, ZN: 8° means 
height at north position at a point 8° latitude upstream from Monterey.) 



Hit R/S | after each entry. 

Enter pressure values for Monterey, Las Vegas, and Eureka. 
Hit 



R/S after each entry. 



10) A data card prompt will occur after all the information has been entered. 
For the first time, use data set #1; the second time, use data set #2, 
etc. 

11) The final number that appears after all data card insertion instructions 
is the probability of average rainfall for the Monterey Peninsula. 

NOTE: If the program ends with "OUT OF RANGE", then the values used are 

either too small or too large. Try another set of values that 

t 

will fit in the stored data sets. 
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D. Examples 

Example A illustrates the use of the program "RAIN" in determining the 
probability of precipitation at Monterey, CA. At 1200Z 25 November 1981 for a 
point 8° latitude upstream from Monterey, the 500 mb heights (in decameters) 
6° latitude to the north, east, south, west and at the center point are 541.0, 
543.0, 562.5, 563.0 and 555.0 respectively. (See Figure 7.) Around Monterey, 
the 500 mb heights are 550.0 dm 6° latitude to the north, 555.0 dm at the 
eastern location, 573.0 dm to the south, 576.0 dm to the west, and 567.0 dm at 
Monterey. These heights are used to obtain the 500 mb relative geostrophic 
vorticity values used in generating forecast variables. The 1200Z sea-level 
pressures at Monterey, Las Vegas and Eureka (1018.8 mb, 1006.2 mb, 1017.3 mb) 
are also required to compute forecast variables. (See Figure 8a.) These 
forecast variables result in a probability of precipitation of 90%. Rain did 
verify during the 24-h period, as noted from the 1200Z 26 November 1981 
observation (Figure 8b.) 
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A 




1200Z 25 November 1981 
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Figure 8a 

Surface observations 
for 1200Z 25 November 1981 



Figure 8b 

Surface observations 
for 1200Z 26 November 1981 



Table II 

Example of the Probability of Precipitation for Monterey, California 
Program. s 



ENTER 


2N : 8 


ENTER- 


ZE:8 


ENTER 


2$: 3 


ENTER 


ZW : 8 


ENTER 


LO • 6 


ENTER 


2H:N 


ENTER 


ZEM 


ENTER 


ZS-.N 


ENTER 


ZW-M 


ENTER 


26 :H 



541. 

543. 

562.5 
563. 
555. 
5*8. 

554. 

568.5 

578.5 

555. 



ENTER PRES MONT 
1.018.8 

ENTER PRES LVGS 
1>006.2800 
ENTER PRES ERKft 
I.- 017..3800 

Kl= 1.018.0000 
X2=-33.0080 
X3= 2.0000 
X4= 14.8080 

INSEPT CRH 13 

Y1 =-20. 0008 

'INSERT CRD 3 

Y2= 25.0000 
80.0000 



Probability of 
Precipitation is 90%. 
Data verified from 
observed case at 
1200Z 26 November 1981. 
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E. Program 



Table III 



Probability 


of 24-h Precipitation for Monterey, 


California 


Program. 


0!*LB‘L “PRIM" 


51 RCL 01 


101 - 


151 GTO OS 


201+LBL 03 


82 “ENTER ZH • 8 ' 


■ 52 t 


182 STO 02 


152 XEQ 00 


292 29.035 


93 PROMPT 


53 CHS 


103 GTO 77 


153 X=Y? 


203 -INSERT CRD 


34 3 TO 33 


54 RCL 30 


184*LBL 37 


154 GTO 09 


204 PVT EH 


35 "ENTER ZE.P- 


55 * 


105 RCL 02 


155 XEQ 03 


205 RDTPX 


33 PROMPT 


56 STO 02 


106 1 


156 X=Y? 


296 FS? PR 


07 RCL 03 


57 -ENTER PRES MONT’ 


107 + 


157 GTO 10 


207 GTO 33 


03 + 


58 PROMPT 


108 STO 82 


158 XEQ 08 


298 FS? 01 


39 3 TO 03 


59 XEQ -ROUND- 


109 GTO 77 


159 X=Y? 


209 GTO 44 


19 "ENTER 23:8- 


68 STO 33 


1104LBL 86 


168 GTO 11 


210 GTO 24 


11 PROMPT 


81 'ENTER PRES LVGS" 


11! RCL 82 


161 XEQ 00 


2!1»LBL 09 


12 RCL 30 


82 PROMPT 


112 -X2=- 


162 X=Y? 


212 20.935 


.13 + 


63 XEQ ■ ROUND - 


113 PCP 


163 GTO 12 


213 -INSERT CRD 


14 5 TO 00 


64 STO 04 


114 sex 


164 XEQ 00 


214 PVIt'H 


15 -ENTER ZW :§■ 


65 -ENTER PRES ERKft* 


115 PDV 


165 X=Y? 


215 RDTPX 


13 PROMPT 


66 PROMPT 


116 RCL 03 


166 GTO 13 


216 FS? 00 


17 RCL 08 


67 XEQ "ROUND- 


117 RCL 01 


167 XEQ 00 


217 GTO 33 


IS + 


68 STO 01 


118 - 


168 X=Y? 


213 FS'** 81 


18 STO 03 


69 RCL 31 


119 XEQ -FIXEVEN- 


169 GTO 14 


219 GTO 44 


23 -ENTER 29 8“ 


70 2 


120 RCL 05 


170 XEQ 00 


228 GTO 24 


21 PROMPT 


71 MOD 


121 STO 00 


171 X=Y? 


22HLBL 10 


22 4 


72 X=8? 


122 *X3=- 


172 GTO 15 


222 20.035 


23 * 


73 GTO 35 


123 PCP 


173 XEQ 00 


223 ’INSERT CRD 


24 CHS 


74 RCL 81 


124 PCX 


174 X=Y? 


224 PV I EH 


25 RCL 30 


75 1 


125 PDV 


175 GTO 16 


225 RDTPX 


23 * 


76 + 


126 RCL 83 


176 FS? 00 


226 FS? 08 


27 STO 03 


77 STO 81 


127 RCL 04 


177 GTO 23 


227 GTO 33 


28 “ENTER ZN:K“ 


78 GTO 35 


123 - 


178 XEQ 88 


228 FS? 01 


28 PROMPT 


79+LBL 85 


129 XEQ ‘FIXEVEN- 


179 X=Y? 


229 GTO 44 


36 STO 31 


39 RCL 01 


138 RCL 05 


130 GTO 17 


230 GTO 24 


31 “ENTER ZE:M- 


31 -XI =■ 


131 STO 04 


181 XEQ 08' 


23!»L6L 11 


32 PROMPT 


32 PCS 


132 -X4=- 


182 X=Y? 


232 20.935 


33 RCL 81 


33 PCX 


133 PCP 


133 GTO IS 


233 “INSERT CRD 


34 + 


34 PDV 


134 PCX 


184 XEQ 80 


234 PVIEW 


35 STO 31 


35 RCL 02 


135 PDV 


135 X=Y? 


235 RDTPX 


36 ‘ENTER ZS:M“ 


36 XEQ -ROUND- 


136 994 


136 GTO 19 


236 FS? 00 


37 PROMPT 


87 STO 02 


137 STO 07 


137 XEQ 00 


237 GTO 33 


33 RCL 01 


88*LBL 77 . 


138 27 


183 X=Y? 


233 FS? 01 


39 + 


89 RCL 62 


139 STO 03 


189 GTO 28 


239 GTO 44 


48 STO 31 


9tf 3 


140 2 


190 XEQ 00 


249 GTO 24 


41 -ENTER ZN:M- 


91 MOD 


141 STO 09 


191 X=Y? 


24i*L8L 12 


42 PROMPT 


92 X=9? 


142 RCL 00 


192 GTO 21 


242 28.035 


43 RCL 31 


93 GTO 86 


143 STO 03 


193 XEQ 00 


243 "INSERT CRD 


44 + 


94 1 


144 RCL 01 


194 X=Y? 


244 PVIEW 


45 STO 31 


95 X=Y? 


145 GTO 40 


195 GTO 22 


245 RDTPX 


46 "ENTER 23 M* 


96 GTO 83 


I46HBL 40 


196 G T Q 23 


246 XROM 23,33 


47 PROMPT 


97 GTO 87 


147 RCL 0? 


197*LBL 00 


247 RCL 01 


48 4 


98*LBL 83 


143 X>Y- 


193 RCL 09 


243 PS? 01 


48 * 


99 RCL 02 


149 GTO 23 


199 + 


249 GTO 44 


50 CHS 


100 1 


150 X=Y? 


230 RTH 


250 GTO 24 
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251+LBl 13 


381 ♦161 18 


351 GTO 51 


401 STO 86 


451 GTO 35 


252 29.835 


382 28.835 


352 XEQ 25 


402 5 


452 RCL 06 


253 “INSERT CRH 6“ 


383 “INSERT CRD If 


353 X=V? 


483 NOB 


453 STO 10 


254 flVIEH 


384 AY I EH 


354 GTO 52 


484 STO 13 


454 “Yl=“ 


255 RDTflX 


305 RDTflX 


355 XEQ 25 


485 X=8’ 


455 PCfl 


258 PV I EM 


386 XRON 28,33 


356 X=Y? 


406 GTO 31 


456 PCX 


257 XRON 28,33 


387 X>8’ 


35? GTO 53 


407 RCL 13 


457 QBV 


258 GTO 33 


388 RCL 12 


358 XEQ 25 


408 1 


458 -2 


259 FS? 81 


389 GTO 24 


359 X=Y? 


489 X=Y’ 


459 STO 87 


269 GTO 44 


318+L8L 19 


368 GTO 54 


418 GTO ?8 


468 18 


261 GTO 24 


311 20.O35 


361 XEQ 25 


411 RCL 13 


461 STO 88 


262»L8L 14 


312 “INSERT CRD !2“ 


362 X=Y? 


412 2 


462 2 


267 28.835 


313 AVIEH 


363 GTO 55 


413 X=Y? 


463 STO 89 


264 “INSERT CRD 7- 


314 RDTflX 


364 XEQ 25 


414 GTO 29 


464 SF 88 


265 AVIEH 


315 GTO 44 


365 X=Y? 


415 RCL 13 


465 RCL 83 


266 RDTflX 


316 XRON 28,33 


366 GTO 56 


416 3 


466 GTO 40 


267 XROn 28,33 


317 GTO 24 


367 XEQ 25 


417 X=Y? 


467 ♦LBL 32 


268 GTO 33 


318+LBL 28 


368 X=Y? 


4 IS GTO 28 


468 RCL 06 


289 FS? 91 


319 20.035 


369 GTO 57 


419 RCL 13 


469 STO 12 


279 GTO 44 


320 “INSERT CRD 13“ 


378 XEQ 25 


428 4 


470 “Y2=‘ 


271 GTO 24 


321 PV I EH 


371 X=Y? 


421 X=Y’ 


471 BCfi 


272*iei_ 15 


322 RDTflX 


372 GTO 58 


422 GTO 27 


472 flCX 


273 28.835 


323 FS? 01 


373 XEQ 25 


427>L6L 27 


473 flDV 


274 “INSERT CRD 3“ 


324 GTO 44 


374 X=Y? 


424 RCL 86 


474 8 


275 AVI EH 


325 CTO 24 


375 GTO 59 


425 l 


475 STO 8? 


276 RDTflX 


326+LBl 21 


376 XEQ 25 


426 ♦ 


476 70 


277 GTO 33 


327 28.035 


377 X=Y’ 


427 STO 06 


477 STO 08 


278 FS’ 33 


328 “INSERT CRD 14* 


378 GTO 60 


428 GTO 31 


473 5 


279*LBL 88 


329 AVIEH 


379 XEQ 25 


429»LBL 28 


479 STO 89 


288 GTO 44 


338 RDTflX 


389 X= v ? 


438 RCL 86 


488 SF 01 


291 GTO 24 


331 AVIEH 


381 GTO 61 


431 2 


481 OF 80 


282»LBL 16 


332 XRON 28.33 


382 XEQ 25 


432 * 


482 RCL 10 


283 28.835 


333 GTO 44 


383 X=Y? 


433 STO 86 


433 GTO 4? 


284 “INSERT CRD 9“ 


334 GTO 24 


334 GTO 62 


434 GTO 31 


484 ♦LBL 23 


285 PVIEH 


335»LBL 22 


335 XEQ 25 


435H.BL 29 


485 FS? 01 


286 RDTflX 


336 28.035 


336 X=Y? 


436 RCL 86 


4:36 GTO “END* 


287 XPOH 28,33 


337 “INSERT CRD 15* 


387 GTO 63 


437 2 


487 “OUT OF RANGE' 


288 RCL 81 


333 AVIEH 


333 XEQ 25 


438 - 


488 AVIEH 


289 FS? 81 


339 RDTflX 


339 X=Y? 


439 STO 86 


489 PTH 


298 GTO 44 


348 GTO d 


398 GTO 64 


448 GTO 31 


4904LBL “END- 


291 GTO 24 


341 XRON 28,33 


391 XEQ 25 


441 ♦LBL 30 


491 FS? 82 


292*L8L 17 


342 RCL 12 


392 X=Y? 


442 RCL 86 


492 GTO 99 


293 29.835 


343 GTO 24 


393 GTO 65 


443 1 


493 SF 02 


294 “INSERT CRD 10* 


344*LBL 24 


394 GTO 26 


444 - 


494 RCL 06 


295 PV I EH 


345 RCL 82 


395*LBL 25 


445 STO 06 


495 1 


296 RDTflX 


346 RCL 88 


396 3 


44S GTO 31 


496 $ T Q 09 


297 XRON 23,33 


347 X=V? 


397 - 


447»LBL 31 


497 XEQ 00 


298»lGl 88 


343 GTO 58 


398 PTH 


448 FS’ 08 


498 SOY’ 


299 GTO 44 


349 XEQ 25 


399HBL 26 


449 GTO 32 


499 GTO 66 


388 GTO 24 


350 X-Y? 


408 XEQ “ROUND- 


458 FS? 81 


580 XEQ 00 
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591 X>Y? 

592 CTO 67 
583 XEQ 89 
594 X>V^ 

585 CTO 98 
598 XEQ 88 
587 X>V? 
598 G'TO 69 
589 XEQ 89 
518 X>Y? 

511 GTO 78 

512 XEQ 80 

513 X>Y? 

514 GTO 71 

515 XEQ 88 

516 X>Y? 

517 GTO 72 

518 XEQ 80 

519 X>Y? 
528 GTO 73 

521 XEQ 88 

522 X>Y? 

523 GTO 74 

524 XEQ 09 

525 X)Y? 

526 GTO 75 

527 XEQ 88 

528 X>Y? 

529 GTO 76 
538 XEQ 08 

531 XEQ 88 

532 X)Y’ 

533 GTO 78 

534 XEQ 80 

535 XEQ 09 

536 X)Y? 

537 GTO 70 

538 XEQ 80 
533 XEQ 89 

540 X>Y? 

541 GTO 30 

542 XEQ 88 

543 XEQ 08 

544 X>Y? 

545 GTO 81 

546 CLX 

547 29 

548 X>Y? 

549 GTO 32 
559 CLX 



551 49 

552 X>Y? 

553 GTO 83 

554 GTO 84 

555 RTH 
556*LBL 99 

557 PCX 

558 ADV 

559 CF 81 
569 CF 92 
561 RTH 
562+LBL 33 



563 


-2 


564 


STO 09 


565 


PCL 84 


566 


PCL 83 


567 


X=Y> 


568 


GTO 58 


569 


XEQ 80 


570 


X = T ' 


571 


GTO 51 


572 


XEQ 80 


573 


X=Y? 


574 


GTO 52 


575 


XEQ 80 


576 


X=Y? 


577 


GTO 53 


578 


XEQ 00 


579 


X=Y> 


538 


GTO 54 


531 


XEQ 00 


532 


X=Y? 


533 


GTO 55 


584 


XEQ 00 


585 


X=Y? 


586 


RCL 26 


537 


GTO 56 


538 


XEQ 80 


539 


X-Y? 


598 


GTO 57 


591 


XEQ 00 


592 


X=Y? 


593 


GTO 58 


594 


XEQ 00 


595 


X=Y? 


596 


GTO 59 


597 


XEQ 00 


593 


X=Y> 


599 


GTO 60 


680 


XEQ 00 



601 


X= Y? 


602 


GTO 61 


683 


XEQ 00 


604 


X=Y? 


605 


GTO 62 


606 


CTO 63 


607 


GTO 26 


688*LBL 35 


609 


PCL 86 


610 


■V jr • 


611 


ACA 


612 


PCX 


613 


PDV 


614 


GTO 23 


615*L8L 44 


616 


-5 


617 


STO 89 


613 


RCL 12 


619 


PCL 98 


629 


X=Y? 


621 


GTO 58 


622 


XEQ 88 


623 


X=Y? 


624 


GTO 51 


625 


XEQ 89 


626 


X=Y? 


627 


GTO 52 


628 


XEQ 80 


629 


X=Y? 


639 


GTO 53 


631 


XEQ 80 


632 


X=Y? 


633 


GTO 54 


634 


XEQ 08 


635 


X=Y? 


636 


GTO 55 


637 


XEQ 89 


633 


X=Y- 


639. 


GTO 56 


648 


XEQ 00 


641 


X=Y? 


642 


GTO 57 


643 


XEQ 00 


644 


X-Y? 


645 


GTO 58 


646 


XEQ 00 


647 


X=Y? 


643 


GTO 59 


649 


XEQ 80 


658 


X=Y'> 



651 GTO 69 

652 XEQ 90 

653 X=x-> 

654 GTO 6l 

655 XEQ 00 

656 X-Y? 

657 GTO 62 

658 XEQ 89 

659 X=Y? 
668 GTO 63 

661 XEQ 08 

662 X=Y? 

663 GTO 64 

664 GTO 65 

665 GTO 26 
666*LBL 59 

667 PCL 28 

668 GTO 26 
669*LBL 51 

678 PCL 21 
671 GTO 26 
672*181 52 

673 PCL 22 

674 GTO 26 
6?5*L8L 53 

676 PCL 23 

677 GTO 26 
678*LBL 54 

679 RCL 24 
638 GTO 26 
681 ♦LBL 55 
682 RCL 25 
633 GTO 26 
684+LBL 56 

635 RCL 26 

636 GTO 26 
687*LBL 57 

638 PCL 27 

639 GTO 26 
690*LBL 58 

691 RCL 23 

692 GTO 26 
693H.BL 59 

694 PCL 29 

695 GTO 26 
696*LBL 60 
697 RCL 30 

693 GTO 26 
699*L6L 61 
700 PCL 31 



701 GTO 26 
702*lBL 62 
703 RCL 32 
794 GTO 26 
705*LBL 83 
706 PCL 33 
797 GTO 26 
7Q8*LBL 64 
789 PCL 34 
719 GTO 26 
711 ♦LBL 65 

712 PCL 35 

713 GTO 26 
714*uBL 66 

715 19 

716 GTO 23 
717*LBL 67 

718 15 

719 GTO 23 
720*LBL 63 
721 25 
72? GTO 23 
723*LBL 69 

724 38 

725 GTO 23 
726 ♦LBL 70 
727 35 
723 GTO 23 
729*LBL 71 

738 40 

73! GTO 23 
732*LBL 72 

733 45 

734 GTO 23 
735*LBL 73 

736 58 

737 GTO 23 
733*LBL 74 

739 55 
748 GTO 23 
74ULBL 75 

742 68 

743 GTO 23 
744*LBL 76 

745 65 

746 GTO 23 
?47*L8L 73 

748 79 

749 GTO 23 
758*LBL 79 



751 75 

752 GTO 27 
753*LBL 80 

754 38 

755 GTO 27 
756 ♦LBL 31 
757 $5 

753 GTO 27 • 

759*L8L 82 

768 98 

761 GTO 23 
762*LBL 83 

763 95 

764 GTO 23 
765* LBL 84 

766 i88 

767 GTO 23 
763*LBL "POUND* 

769 FIX 8 
778 PhD 

771 FIX 4 

772 RTN 

773*LBL "FIXEVEH" 

774 STO 85 

775 2 

776 HOD 

777 X=07 

778 RTN 

779 PCL 85 
738 1 

731 + 

732 870 85 
783 PIN 
734 .END. 
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IV. Contrail Formation 



A. Objective 

To determine whether or not contrails will form at a certain level in the 
atmosphere based on the curves from the Contrail Forecasting Manual (Chief of 
Naval Operations, 1964). 

B. Principles 

The critical temperature determines if contrails will form at a level. 
A least squares approximation is given for the critical temperature: 

Ter it = ai + aa In p + a3 (In p) 2 + a4 RH + as (RH) 2 ; 

( 

where: 

p is the pressure level; 

RH is the relative humidity; 

ax = -90.4994; as = 3.4232; a 3 = 0.5587; a 4 = -0.0372; as = 0.0012. 

\ 

This equation was determined from the curves in Figure 9 of the Contrails 
Forecasting Manual (Chief of Naval Operations, 1964). 

If the ambient temperature is less than the critical temperature, then 
contrails will form. If the ambient temperature is greater than the critical 
temperature, then contrails will not form. 

There is a narrow range of temperatures where the relative humidity is 
important. If the dew point is known, the relative humidity can be calculated 
by using the definitions of , saturation vapor pressure (es) and relative 
humidity (R.H.). 

es(T) = 6.11 exp (L/Rw (1/273 - 1/T)) 

R.H. = 100 * (e/ea) 

where es (T) is the saturation vapor pressure 
e is the actual vapor pressure 
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T is the temperature (°K) 

L is the latent heat of evaporation (J/kg) 

Rw is the specific gas constant for pure water vapor (J/kg °K) 

If the dew-point temperature is unknown, the relative humidity is 
estimated and a +2°C error margin is assumed when computing the critical 
temperature. If the ambient temperature is within +2°C of the critical 

temperature, a "probably" will precede the "contrails" or "no contrails" 
message. The following criteria are used to estimate the relative humidity: 

Tamb < Ter it — > (RH = 0) — > contrails 
Tamb < Ter it — > (RH = 100) — > no contrails 
p < 225 — > (stratosphere) — > RH — 0 % 
p > 300 — > (not in upper troposphere) — > RH — 40% 

225 < p < 300 — > (upper troposphere) — > RH — 40 % UNLESS 
cirrus are at this level — > RH — 60 % 
flow is known and from 

moist region — > RH — 60% 
dry region — > RH — 0 % 
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c. 


User’s 


Guide 






1) 


Make sure both the calculator and 


printer are on and 




set to 


the MAN mode. 






2) 


XEQ 


alpha SIZE alpha 






3) 


Prompts SIZE ; Type 100. 






4) 


Load all three cards, sides 1-6. 






5) 


After ’ 


’WORKING" signal, program is 


ready for 


execution 


6) 


XEQ 


alpha TRAILS alpha 






7) 


Program asks about CHANGING FLAGS. 


For the 


first run 



8 ) 

9) 



flags as they are and continue by hitting the | R/S key. 
Enter a pressure level in mb. 

Press 



R/S key. (Hit jR/S| after each entry.) 



10) Enter a temperature in °C. 

11) A prompt for a dew-point temperature will be displayed. If Td is 
known, type a 1; if Td is unknown, type a 0. 

12) The program will then decide whether or not contrails can be formed at 
that particular level depending upon the criteria discussed on the 
preceding pages. If the dew-point temperature is unknown, there might be 
more "yes" or "no" questions depending upon the pressure level and 
ambient temperature. Answer "l” for "yes" and "0" for "no.” 

13) The program will determine whether or not contrails 1 will form at a 
certain pre-established pressure level. 

14) After the answer is displayed, execution can continue for another level 
by pressing | R/S ] 
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By changing the flags, two options are possible. Setting Flag 1 will 
make the prompts shorter and less specific. Setting Flag 2 means that 
the program assumes that the dew-point temperature is known and skips the 



questions concerning relative humidity estimation. 

To change a flag: a) Press gold button and hei key. 

b) Program prompts SF . 

c) Enter either (or both): 01 - short prompt 

02 - known Td 

d) To clear flags, press gold button and HE key. 

e) Program prompts CF . 

f) Enter 01 or 02 to clear any flag that was set. 
The program assumes all flags are clear each time it is first executed. 



D. Examples 



Examples A-E in 


Table IV 


illustrate 


the 


use of the program "TRAILS" 


in 


deciding whether 


contrails 


will 


form 


from 


upper-level sounding data. 


In 


example A, where 


no flags 


have 


been 


set, 


a pressure of 450 mb 


and 



temperature of -35°C are entered in the "TRAILS'' program. The dew-point is 
unknown in this case so a "0" is entered. The program decides that contrails 
will not be possible with these conditions which agrees with the contrail 
curve diagram. (See Figure 9.) Example B is similar to A (no flags set and 
dew point unknown) , but contrails will probably form with the particular 
upper-level sounding data. In example C, no flags are set, but the dew point 
is known. With a pressure of 250 mb, a temperature of -60.0°C, and dew point 
of -62.0°C, contrails will form. Examples D and E have flags set, which 
result in shortened prompts for both cases, and for example E, the dew point 
is assumed to be known. Both examples D and E illustrate that contrails 
will not form in these environmental conditions. 
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Table IV 

Examples of the Contrail Formation Program. 



A 




B 








P=' 7 HB 




P=? MB 




250.0000 


*** 


450.0000 


*** 


T=? C 




t='> r 




-54.0000 


*** 


-35.0000 


*** 


DO YQU KNOW 




DO YOU KNOW 




DEWPOINT, TD? 




DEWPOINT. W 




YES=ONE 




YES=ONE 




N0=2EP0 




NO=ZERG 




KNOW TD' 7 




KNOW TD? 




0.0000 


*** 


0.0000 


*** 


ARE THERE 




NO CONTRAILS 




CIRRUS AT 
THIS LEVEL' 7 




No flags set 
unknown 




VES=ONE 

NO=ZERO 

CIRRUS' 7 








1.0000 


*** 






P R0BA8LY 

CONTRAILS 





P =' 7 .MB 

2 25.8809 

!=■? c 

- 68.0009 
DO you KNOW 
DEWPOINT. TD 7 
Yc$=0NE 
N0=2ER0 
KNOW ID? 

1,0000 

TD=? C 

- 62.8808 

CONTRAILS 



*** 



■***. 



*** 

*** 



No flags set, Tj unknown 

D 



P=? HP. 

T=? C 
KNOW ID? 
CIRRUS? 

KNOW FLOW' 7 
MST OR DRY? 

1 

NO CONTRAILS 



300.0000 
-40.0000 
8.0000 
0.0000 
1.0000 
. 8080 



*** 

*** 

*** 

**■* 

*** 

*** 



No flags set 
T ^ known 



E 



P=? MB 

050. 0000 

T= 7 C 

-18.0000 



TD=? C 

- 20.0000 
NO CONTRAILS 



Flag 01 set 
Unknown T d 



*** 

*** 

**■* 



Flag 02 set 
T ^ known 
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Figure 9 

Illustration of contrail decisions: curves taken from the Chief of 
Naval Operations C1964). 
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E. Program 





Table V 


Contrail Formation 


Program. 




9! ♦Let "TRAILS' 


51 GTO 11 


101 STO 82 


151 XROH 28,33 


20! T ONt 3 


92 CF 9l 


52 -BO YOU KNOW** 


102 FS? 81 


152 X=0? 


202 "TD=? C- 


93 CF 82 


53 AVIEH 


183 GTO 12 


153 GTO A 


203 FS’ 21 


94 -CHANGE FLAGS'*' 54 PSE 


184 -ARE THERE- 


154 FS? 81 


204 PSA 


05 PROMPT 


55 -DEHPOIHT, TP 


?■ 195 OVIEH 


155 GTO 14 


285 FS? 2! 


96+LSL 19 


56 AVIEH 


186 PSE 


156 -FROH HOIST- 


206 PRX 


07 TONE l 


57 PSE 


107 -CIRRUS AT- 


157 AVIEH 


207 ‘XROH 28,33 


08 *P=? HB- 


58 -YESsOHE- 


198 AVIEH 


158 PSE 


288 273 


OS FS? 21 


59 AVIEH 


189 PSE 


159 -REGION- 


289 * 


10 PRA 


68 PSE 


110 -THIS LEVEL?" 


168 AVIEH 


210 XES ■£$* 


11 PROMPT 


61 -NO=ZERO- 


Ill AVIEH 


161 PSE 


211 RCL 02 


12 FS? 21 


62 AVIEH 


112 PSE 


162 "= ONE" 


212 / 


13 F'RX 


63 PSE 


113 "YES=ONE“ 


163 AVIEH 


213 109 


14 STO 00 


64*LBL 11 


114 AVIEH 


164 PSE 


214 * 


15 LN 


65 -KNOW TD?- 


115 PSE 


165 -FROM PRY- 


215 STO 02 


16 STO 01 


66 FS? 21 


116 -HO=ZERO- 


166 AVIEH 


216 RTH 


17 TONE 2 


67 PRA 


117 AVIEH 


167 PSE 


217+LBL -ES- 


18 “T=> IS- 


63 XROH 28,21 


118 PSE 


163 •= ZERO’ 


218 1/X 


IS FS? 21 


69 XROH 28,33 


li9*LBL 12 


169 AVIEH 


219 CHS 


'28 PRA 


TO CHS 


120 -CIRRUS’- 


178 PSE 


220 273 


21 PROMPT 


71 1 


12! FS? 21 


171 ♦LBL 14 


221 1/X 


22 FS’ 21 


72 X=Y? 


122 PRA 


172 -HST OR PRY’- 


222 ♦ 


23 PRX 


73 GTO 01 


123 FS? 21 


173 FS? 21 


223 2580 


24 STO 04 


74 0 


124 XROH 23,33 


174 PRA 


224 * 


25 FC? 02 


75 STO 82 


125 PRX 


175 XROH 28,33 


225 0.461 


28 GTO 03 


76 XES -TC- 


126 X=Y? 


176 5 


226 / 


27*L8L 01 


77 RCL 84 


127 GTO A 


177 PRX 


227 EtX 


28 XEQ -RH- 


78 RCL 83 


128 XROH 23,33 


178 * 


228 6 


2? XEQ “TO- 


79 X>Y? 


129 48 


179 GTu 39 


229 * 


30*LBL 02 


80 GTO C 


138 STO 82 


180 RCL 81 


230 RTN 


31 PCL 04 


81 188 


131 FS? 01 


181 X<0? 


231 ♦LBL -TC- 


32 RCL 93 


82 STO 92 


132 GTQ 12 


182*LBL ? 


232 -90.4994 


33 X>Y? 


83 XEQ -TC* 


133 -DO YOU- 


183 XEQ -TC- 


233 STO 03 


34 GTO C 


84 RCL 83 


134 AVIEH 


184 pci 83 


234 .5587 


35+LBL B 


85 RCL 04 


135 PSE 


185 RCL 04 


235 RCL 01 


36 *NO CONTRAILS' 


86 X>Y’ 


136 -KNOW FLOH?- 


186 - 


236 * 


37 TONE 5 


87 GTO B 


137 AVIEH 


187 ABS 


237 3.4232 


38 AVIEH 


88 0 


138 PSE 


133 2 


238 ♦ 


39 ADV 


39 STO 82 


139 ‘YESsOHE- 


189 X<=Y? 


239 RCL 91 


49 STOP 


98 RCL 00 


140 AVIEH 


198 GTO 82 


240 * 


41 GTO 10 


91 225 


141 PSE 


191 -PROBABLY- 


241 ST+ 03 


42*LBL C 


92 X)Y? 


142 -NO=ZERO- 


192 AVIEH 


242 .0012 


43 -CONTRAILS- 


93 GTO A 


143 AVIEH 


193 PSE 


243 RCL 02 


44 TONE 8 


94 48 


144 PSE 


194 GTO 02 


244 


45 AVIEH 


95 STO 82 


145+LBL 13 


195*lBL -RH* 


245 -.0372 


46 ABV 


96 309 


146 ‘KNOW FLOH?' 


196 RCL 04 


246 ♦ 


47 STOP 


97 RCL 09 


147 FS? 21 


197 273 


247 RCL 02 


48 GTO 10 


98 X>Y? 


148 PRA 


198 + 


248 * 


49»LBL 93 


99 GTO A 


149 FS? 21 


199 XEQ -ES- 


249 ST+ 93 


59 FS? 01 


100 60 


158 PRX 


208 STO 02 


250 RTN 

251 .END. 
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V. Determination of Lifting Condensation Level 

A. Objective 

Given a surface temperature, surface pressure and specific humidity, this 
program calculates the lifting condensation level. 

B. Principles 

The lifting condensation level (LCL) is the level to which a parcel of 
air can be lifted dry adiabatically before it becomes saturated. During the 
lifting process, the potential temperature of the air parcel and the 
saturation mixing ratio remain constant. The actual mixing ratio decreases 
and eventually equals the saturation mixing ratio. The lifting condensation 
level is obtained when the actual vapor pressure equals the saturation vapor 
pressure (e = es ) . Integrating the Clausius-Clapeyron equation from 273 K to 
air temperature (T) 




des 

e s 




dT 

T 2 



where es is the saturation vapor pressure; L is the latent heat of 
condensation (2.5 x 10 6 J/kg) ; Rw is the gas constant for water vapor (462 
J/kg °K); and T is the temperature (°K), results in the following expression 
for the saturation vapor pressure: 



es =6.11 exp (-— (~ - h). 

Rw 273 T 

Combining the hydrostatic equation and the ideal gas law and expressing the 
temperature as To - Y z results in the following expression relating pressure 
(P) to height (z): 

P . p„ 

To 

I 
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where y is the dry adiabatic lapse rate, To is the standard atmospheric 
temperature, w is the specific humidity and Po is the standard atmospheric 
pressure, (w = .622 (e/P) ==> P = .622 e/w) 

An expression for the actual vapor pressure (e) is obtained by using the 
definition of specific humidity and substituting the above expression in for 
P. The lifting condensation level is obtained when the following two 
equations are equal to each other: 



After a series of substitutions, using Taylor series expansion of 
In x = x-1 and several algebraic manipulations, the following equation is 
derived giving the lifting condensation level: 



_wPo (To - yz ) K 
. 622 To 





zlcl = To (102.041 - 14.6429 (23.3058 - ln( 



wPo 

3. 80042 } 



[ln( 3.80042 ) 



- 23.3058 2 



75418.2 

To 



] 2 )) 



Only surface data values of temperature (To), pressure (Po) and specific 
humidity (w) are required to find the LCL. 



c. 



1 ) 

2 ) 

3) 

4) 

5) 

6) 

7) 

8 ) 

9) 

10 ) 
11 ) 



User’s Guide 

Turn on calculator and set printer to the NOR mode. 



XEQ 



LCL 



Prompt reads, SFC TEMP °C? 
Type in a surface temperature. 
Press 



R/S to enter the temperature. 



Prompt is SFC PRES MB? 
Enter a surface pressure. 



R/S 

Prompt reads, SP HUMID G/G? 
Enter a specific humidity value. 
R/S 



The lifting condensation level will appear with its value in meters. If 
the LCL occurs at the surface, the word SURFACE will be displayed. 



D. Examples 

Examples A, B and C in Table VI illustrate the use of the program "LCL" 
in estimating the lifting condensation level from a surface temperature, 
pressure and specific humidity. In example A, a surface temperature of 9°C, 
a surface pressure of 1010 mb, and a specific humidity of .006 g/g are entered 
in the "LCL" program. The data result in a lifting condensation level of 
318.9m which is verified by a thermodynamic diagram. (See Figure 10.) In 
example B, the lifting condensation level occurs at the surface and in example 
C, the LCL is at 637.4 m. 
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Table VI 

Examples of the Lifting Condensation Level Program. 



A 



B 



S*C TEH? C? 

9,0000 

S C C PRES MB' 7 

1 / 010.0000 
SP HUH ID G/G? 

. 0060 

l CL/ H/=31095‘? . 



SFC TEMP C' 7 

. 5000 

SFC PRES MB? 

1 / 000.0000 
SP HUH ID G/G’ 

.0040 

SURFACE 




SFC TEMP C? 

6.0000 

SFC PRES MB? 

1 / 020.0000 
SP HUH ID G/G? 

.9040 

LCL, H- =637.. 3732 
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Figure 10 

Lifting condensation levels determined by surface temperature, pressure 
and specific humidity. 
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E. Program 



Table VII 

Lifting Condensation Level Program. 



01+L8L "LCl" 

02 23.3058 
0? '2 Tfi 25 

04 -SPG TEMP 

05 PROMPT 

88 273.15 
87 + 

98 8 TO 21 

89 “SFC PRES NS ?" 
18 PROMPT 

11 STO 22 

12 *SP M1JNJ3 Q 
17 PROMPT 

14 STO 23 

15 RCL 22 

16 * 

17 3.38042 
IS / 

1? in 
28 STO 24 

21 RCL 23 

22 - 

23 >; + 2 

24 75418.2 

25 RCL 21 

26 / 

27 - 

23 SORT 



29 STO 26 

30 CPS 

31 RCL 24 

32 CHS 

33 * 

34 RCL 25 

35 * 

36 14.64V? 

37 t 

38 182.041 

39 XOY 

40 - 

4! RCL 21 

42 * 

43 STO 38 

44 A 

45 XOY 

46 X ! "=Y? 

47 CTO 'SRC 
43 "Lit. H..=’ 

49 fiRCL 38 

50 hVIEH 

51 CTO “END* 
52+LEL *S C C* 

53 -SOPPflCE’ 

54 flVIEH 
55+LBL “EH?* 
56 .END. 
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